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Efficient methods and apparatus for resampling three-dimensional datasets 



(57) In one form, the present invention is a method 
for rapidly re-sampling a three-dimensional volume of 
data (52) to permit efficient interactive viewing and visu- 
alization of the data is provided in which a set of parallel 
line segments (56, 58, 60, 62, 64, 66, 68) defining a sur- 
face (50) having a constant slope in at least one direc- 
tion (A) and intersecting the three-dimensional volume 
are selected and the data is processed at points (132, 
140, 146, 150, 152, 154, 156, 158, 160, 162, 166) on 
the parallel line segments to generate a reformatted 
image of the data. A corresponding apparatus embodi- 
ment is provided. Processing of data can proceed 
quickly because it is not necessary to check the parallel 
line segments during processing of the points (158, 
160, 162, 164, 166) on the segments to determine 
whether the points are within a given data volume. 
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Description 

[0001] This invention relates generally to methods 
and apparatus for re-sampling of data and more partic- 
ularly to methods and apparatus for re-sampling three 5 
dimensional datasets for visual processing and display 
of medical data. 

[0002] In at least one known computed tomography 
(CT) imaging system configuration, an x-ray source 
projects a fan-shaped beam which is collimated to lie w 
within an X-Y plane of a Cartesian coordinate system 
and generally referred to as the "imaging plane". The x- 
ray beam passes through the object being imaged, such 
as a patient. The beam, after being attenuated by the 
object, impinges upon an array of radiation detectors, is 
The intensity the attenuated beam radiation received 
at the detector array is dependent upon the attenuation 
of the x-ray beam by the object. Each detector element 
of the array produces a separate electrical signal that is 
a measurement of the beam attenuation at the detector 20 
location. The attenuation measurements from all the 
detectors are acquired separately to produce a trans- 
mission profile. 

[0003] In known third generation CT systems, the x- 
ray source and the detector array are rotated with a gan- 25 
try within the imaging plane and around the object to be 
imaged so that the angle at which the x-ray beam inter- 
sects the object constantly changes. A group of x-ray 
attenuation measurements, i.e., projection data, from 
the detector array at one gantry angle is referred to as a 30 
"view*. A "scan" of the object comprises a set of views 
made at different gantry angles, or view angles, during 
one revolution of the x-ray source and detector. In an 
axial scan, the projection data is processed to construct 
an image that corresponds to a two dimensional slice 35 
taken through the object. One method for reconstructing 
an image from a set of projection data is referred to in 
the art as the filtered back projection technique. This 
process converts the attenuation measurements from a 
scan into integers called "CT numbers" or "Hounsfield 40 
units", which are used to control the brightness of a cor- 
responding pixel on a cathode ray tube display. These 
integers form a volume of medicaJ data. 
[0004] One popular visualization technique for 
medical applications is to resample a volume of medical 45 
data acquired from a CT or MRI imaging system along 
a plane of arbitrary slope. This visualization technique is 
known as a reformat. Volume re-sampling along a plane 
is commonly called a reformat of a given volume. Refor- 
matting allows viewing of cross-sections of pathology so 
and provides the ability to take accurate measurements 
of strictures, aneurysms, areas, and distances, among 
other things. 

[0005] Determining the placement of a re-sampling 
plane could be facilitated by interactive placement. 55 
Interactive placement would provide for the display of a 
first reformat plane image followed by the production 
and display of a new reformat plane image, where the 



new reformat plane is incrementally different in scale, 
rotation, or translation from the first reformat plane. 
Interactive placement, however, demands good re-sam- 
pling performance of the data, usually requiring about 
10 reformat images per second to be generated. With- 
out adequate performance, interactive placement is dif- 
ficult to control. Systems known to date have not been 
capable of providing adequate re-sampling perform- 
ance. 

[0006] Accordingly, it would be desirable to rapidly 
and efficiently resample a three-dime nsionai volume of 
data to permit efficient interactive viewing and visualiza- 
tion of medical data. It would also be desirable to pro- 
vide efficient clipping of re-sampling surfaces. In 
addition, it would be desirable to provide efficient man- 
agement of optimized re-sampling functions, where 
large sets of highly optimized and specialized functions 
are used for re-sampling. 

[0007] According to a first aspect of the invention, 
there is provided a method for rapidly re-sampling a 
three-dimensional volume of data to permit efficient 
interactive viewing and visualization of the data, the 
method comprising: selecting a set of parallel line seg- 
ments defining a surface having a constant slope in at 
least one direction and intersecting the three-dimen- 
sional volume; and processing the data at points on the 
parallel line segments to generate a reformatted image 
of the data. 

[0008] Selecting a set of parallel line segments may 
comprise the step of clipping ends of the line segments 
so that the clipped line segments lie wholly within the 
three-dimensional volume, and the processing step 
comprises processing the data along at points on the 
clipped line segments. 

[0009] Selecting a set of parallel line segments may 
comprise the step of selecting a set of coplanar parallel 
line segments or may comprise the step of selecting a 
set of parallel line segments in a surface having a pair of 
orthogonal axes, the surface having curvature along 
exactly one of the pair of orthogonal axes. 
[0010] Clipping ends of the line segments may 
comprise the steps of: determining a first set of faces of 
the volume for clipping a first endpoint of each of the line 
segments and a second set of faces of the volume for 
clipping a second, opposite endpoint of each of the line 
segments; and for each line segment, analyzing the first 
endpoint of the line segment and the first set of faces to 
clip at a first clipping point, and analyzing the second 
endpoint of the line segment and the second set of 
faces to clip at a second clipping point 
[001 1] Determining a first set of faces and a second 
set of faces may comprise the step of determining a first 
set of three faces and a second set of three faces. 
[0012] Processing of the data may comprise the 
steps of selecting a specialized sampling function corre- 
sponding to a slope of the parallel line segments and 
computing the selected specialized sampling of the data 
at points along the clipped line segments. 
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[0013] Selecting a specialized sampling function 
corresponding to a slope of the parallel line segments 
may comprise the step of selecting one of a plurality of 
specialized sampling functions in accordance with a 
slope of the parallel line segments. 
[0014] Selecting one of a plurality of specialized 
sampling functions may comprise the step of selecting a 
specialized sampling function having a selected interpo- 
lation method from among several specialized sampling 
functions corresponding to a slope of the parallel line 
segments and that provide different interpolation meth- 
ods. 

[0015] Selecting one of a plurality of specialized 
sampling functions may comprise the step of selecting a 
specialized sampling function having a selected voxel 
addressing scheme from among several specialized 
sampling functions corresponding to a slope of the par- 
allel line segments and that provide different voxel 
addressing schemes. 

[0016] The method may further comprise the step 
of generating the sampling functions from a template. 
[0017] Automatically generating the sampling func- 
tions may comprise the step of generating the sampling 
functions using class-structured sets of macros. 
[0018] Selecting a specialized sampling function 
having a selected interpolation method may comprise 
the step of selecting a trilinear interpolation, nearest 
neighbor, or process sampling function. 
[0019] Computing the selected specialized sam- 
pling of the data at points along the clipped line seg- 
ments may comprise the step of applying a step-and- 
interpolate iteration to the sample data, and the applied 
step-and-interpolate iteration may be free of bounds 
checking. 

[0020] The method may further comprise the step 
of selecting a baseline segment forming one edge of a 
sampling region. 

[0021] According to a second aspect of the Inven- 
tion, there is provided a system for rapidly re-sampling a 
three-dimensional volume of data for efficient interactive 
viewing and visualization of the data, said system con- 
figured to: select a set of parallel line segments defining 
a surface having a constant slope in at least one direc- 
tion and intersecting the three-dimensional volume; and 
process the data at points on the parallel line segments 
to generate a reformatted view of the data. 
[0022] The system may be further configured to clip 
ends of the parallel line segments so that the clipped 
line segments lie wholly within the three-dimensional 
volume, and to process the data at points on the clipped 
line segments. 

[0023] The system may be further configured to 
select a set of parallel line segments in a surface having 
a pair of orthogonal axes, the surface having curvature 
along exactly one of the pair of orthogonal axes. 
[0024] The system may be further configured to 
select a specialized sampling function corresponding to 
a slope of the parallel line segments and to compute the 



selected specialized sampling of the data at points 
along the clipped line segments. 
[0025] The system may be further configured to 
determine a first set of faces of the volume for clipping a 

5 first endpoint of each of the line segments and a second 
set of faces of the volume for clipping a second, oppo- 
site endpoint of the line segments, to analyze, for each 
line segment, the first endpoint of the line segment and 
the first set of faces to clip at a first clipping point, and to 

io analyze the second endpoint of the line segment and 
the second set of faces to clip at a second clipping point. 
[0026] The system may be further configured to 
determine a first set of three faces as the first set of 
faces of the volume, and to determine a second set of 

15 three faces as the second set of faces of the volume. 
[0027] The system may be further configured to 
select a set of coplanar parallel line segments. 
[0028] The system may be further configured to 
select one of plurality of specialized sampling functions 

20 in accordance with a slope of the parallel line segments. 
[0029] The system may be further configured to 
select a specialized sampling function having a selected 
interpolation method from among several specialized 
sampling functions corresponding to a slope of the par- 

25 allel line segments and that provide different interpola- 
tion methods. 

[0030] The system may be further configured to 
select one of a plurality of specialized sampling func- 
tions having a selected voxel addressing scheme from 

30 among several specialized sampling functions corre- 
sponding to a slope of the parallel line segments and 
that provide different voxel addressing schemes. 
[0031] The system may be further configured to 
generate sampling functions from a sampling function 

35 template. 

[0032] The system may be further configured to 
generate sampling functions using class-structured sets 
of macros. 

[0033] The system may be further configured to 
40 select a specialized sampling function having a selected 
interpolation method from a group of sampling methods 
including trilinear interpolation, nearest neighbor, and 
process sampling functions. 

[0034] The system may be further configured to 
45 apply a step-and-interpolate iteration to the sample data 
as part of the computation of the selected specialized 
sampling of the data at points along the clipped line seg- 
ments. 

[0035] The system may be further configured so 
so that the step-and-interpolate iteration is free of bounds 
checking. 

[0036] The system may be further configured to 
select a baseline segment forming one edge of a sam- 
pling region. 

55 [0037] In one embodiment of the invention a 
method for rapidly re-sampling a three-dimensional vol- 
ume of data to permit efficient interactive viewing and 
visualization of the data in which a set of parallel line 
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segments defining a surface having a constant slope in 
at least one direction and intersecting the three-dimen- 
sional volume are selected and the data is processed at 
points on the parallel line segments to generate a refor- 
matted image of the data. A corresponding apparatus 
embodiment is provided. Using the methods and appa- 
ratus embodiments of this invention allows processing 
of data to proceed quickly because it is not necessary to 
check the parallel line segments during processing of 
the points on the segments to determine whether the 
points are within a given data volume. Furthermore, the 
methods and apparatus embodiments allow rapid and 
efficient re-sampling of a three-dimensional volume of 
data, for example, medical data, to be achieved through 
use of sets of highly optimized and specialized re-sam- 
pling functions. 

[0038] The invention will now be described in 
greater detail, by way of example, with reference to the 
drawings, in which:- 

Figure 1 is a pictorial view of a CT imaging system. 

Figure 2 is a block schematic diagram of the system 
illustrated in Figure 1. 

Figure 3 is a drawing illustrating an intersection of a 
re-sampling plane with a data volume represented 
by a cube. 

Figure 4 is a drawing illustrating the intersection of 
a more general curved surface with a data volume 
represented by a cube. 

Figure 5 is a flow chart of an embodiment of a line 
segment clipping procedure. 

Rgure 6 is a drawing illustrating a simplified, two- 
dimensional example of an embodiment of a line 
segment clipping procedure as applied to a line 
segment intersecting a square, two-dimensional 
data volume. 

Rgure 7 is a schematic representation of an 
embodiment of a decision tree useful for selecting a 
specialized sampler. 

Rgure 8 is a drawing of a simplified illustration of 
planar re-sampling in which a baseline segment 
forms one edge of a sampling region and line seg- 
ments normal to the baseline segment provide 
scanline segments. 

[0039] Referring to Rgures 1 and 2, a computed 
tomograph (CT) imaging system 1 0 is shown as includ- 
ing a gantry 1 2 representative of a "third generation" CT 
scanner. Gantry 1 2 has an x-ray source 1 4 that projects 
a beam of x-rays 16 toward a detector array 18 on the 
opposite side of gantry 12. Detector array 18 is formed 



by detector elements 20 which together sense the pro- 
jected x-rays that pass through a medical patient 22. 
Detector array 18 may be fabricated in a single slice or 
multi-slice configuration. Each detector element 20 pro- 

s duces an electrical signal that represents the intensity 
of an impinging x-ray beam and hence the attenuation 
of the beam as it passes through patient 22. During a 
scan to acquire x-ray projection data, gantry 12 and the 
components mounted thereon rotate about a center of 

70 rotation 24. 

[0040] Rotation of gantry 12 and the operation of x- 
ray source 1 4 are governed by a control mechanism 26 
of CT system 1 0. Control mechanism 26 includes an x- 
ray controller 28 that provides power and timing signals 

15 to x-ray source 1 4 and a gantry motor controller 30 that 
controls the rotational speed and position of gantry 12. 
A data acquisition system (DAS) 32 in control mecha- 
nism 26 samples analog data from detector elements 
20 and converts the data to digital signals for subse- 

20 quent processing. An image reconstructor 34 receives 
sampled and digitized x-ray data from DAS 32 and per- 
forms high speed image reconstruction. The recon- 
structed image is applied as an input to a computer 36 
which stores the image in a mass storage device 38. 

25 [0041] Computer 36 also receives commands and 
scanning parameters from an operator via console 40 
that has a keyboard. An associated cathode ray tube 
display 42 allows the operator to observe the recon- 
structed image and other data from computer 36. The 

30 operator supplied commands and parameters are used 
by computer 36 to provide control signals and informa- 
tion to DAS 32, x-ray controller 28 and gantry motor 
controller 30. In addition, computer 36 operates a table 
motor controller 44 which controls a motorized table 46 

35 to position patient 22 in gantry 12. Particularly, table 46 
moves portions of patient 22 through gantry opening 48. 
[0042] A scan of patient 22 results in the collection 
of a volume of medical data. Interactive volume re-sam- 
pling of this data allows viewing cross-sections of 

40 pathology to take accurate measurements of strictures, 
aneurysms, areas, and distances. In one embodiment, 
the invention provides a fast and efficient clipping 
method that is applicable to re-sampling surfaces that 
can be modeled as a set of line segments of identical 

45 slope. These surfaces include planes and other sur- 
faces having two orthogonal axes and curvature along 
one of the axes, but not the other. One example of such 
a surface is a Z-plane, with a displacement along the x 
axis, or along the y axis. These types of surfaces can 

so have key clipping aspects determined for all line seg- 
ments on the surface, accelerating the clipping process 
for each line segment. 

[0043] More generally, the intersection of a re-sam- 
pling plane 50 with a data volume 52, represented, for 
55 example, by a cube, is illustrated in Rgure 3. Dashed 
lines (e.g., 54) are shown to represent lines on which 
sampling points interior to volume 52 lie. In one embod- 
iment, re-sampling at points in data volume 52 begins 
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by clipping each of the line segments in plane 50 until 
sampling points are interior to the volume 52, or a line 
segment being analyzed is found to be entirely outside 
volume 52. Line segments 56, 58, 60, 62, 64, 66, and 68 
lying within plane 50 are representative of the types of 
line segments to be analyzed for clipping. Some of 
these line segments (56, 58, 66, 68) lie entirely outside 
the volume. Others (60, 62, 64) lie partially within the 
volume and are clipped for re-sampling. The portion of 
data volume 52 along which data points are re-sampled 
lie within a sheet 64. In Figure 3, parallel line segments 
56, 58, 60, 62, 64, 66, and 68 of equal slope define an 
axis A of plane 50. As a result of the re-sampling sur- 
face being planar, there is also no curvature along 
orthogonal axis B of plane 50. 

[0044] In another embodiment, a curved re-sam- 
pling surface is examined. For example, figure 4 repre- 
sents the intersection of a curved re-sampling surface 
72 with data volume 52. Re-sampling surface 72 is 
curved, but is still characterized or defined by parallel 
line segments 74, 76, 78, 80, 82, 84, and 86 having 
equal slope. The slope of these parallel line segments 
defines an axis A of re-sampling surface 72. There is no 
curvature of re-sampling surface 72 in directions paral- 
lel to axis A, but curvature is allowed, and is, in fact, 
present, in the direction of an axis B, which is orthogo- 
nal to axis A. In the example shown in Figure 4, the re- 
sampling surface is a curved sheet 88 lying entirely 
within volume 52. 

[0045] Segment ends are clipped back until re-sam- 
pling points are interior to volume 52 or the segment is 
found to be entirely outside the volume. In one embodi- 
ment, the process is performed in two iterates, one for 
each end of the line segment. The process starts with 
one extreme sample point along the line segment and 
tests each sample point until either the other extreme 
end of the line segment is encountered or a sample 
point interior to the volume is found. Because the seg- 
ment endpoints are clipped in the beginning, step-and- 
interpolate iterations that process the sample points 
within the data volume can be performed free of bounds 
checking. In other words, it is not necessary to check 
each point within the step-and-interpolate iterations to 
determine whether the point being processed is inside 
or outside the data volume. The efficiency of the refor- 
matting operation is increased because this check is 
eliminated from the iterations. 

[0046] Figure 5 is a flow chart of a clipping proce- 
dure that is used in one embodiment. Execution of the 
procedure starts at block 90, where one of the two 
opposite endpoints of a line segment of the sampling 
plane is selected. At block 92, a test is performed to 
determine whether the selected point is within a bound- 
ary of data volume 52 represented by the lower faces of 
the volume. This boundary corresponds, for example, to 
bottom face 94, right face 96, and front face 98 of vol- 
ume 52 in Figures 3 and 4, but different faces appropri- 
ate for this test are used, depending upon an orientation 



of a sampling surface, for example, 50 or 72. If the point 
is not within the lower faces, i.e., within the data volume, 
execution proceeds to block 100. Otherwise, execution 
continues at block 102. 

5 [0047] Block 1 00 is executed only if the re-sampling 
point under consideration is outside data volume 52. In 
this case, the output point, i.e., the visual representation 
point or pixel corresponding to the re-sampling point, is 
filled with "background" data to indicate that the output 

w point is outside data volume 52 and thus, no medical 
data is available. Execution then proceeds to block 1 04, 
in which a test is performed to determine whether the 
opposite end of the line segment has been reached. If 
so, it is not necessary to check the line segment any f ur- 

15 ther, because the line segment lies entirely outside of 
data volume 52. Execution thus returns to the calling 
procedure at block 106, which either checks another 
line segment of the re-sampling surface or proceeds to 
other processing. Otherwise, a step is taken along the 

20 current line segment towards the opposite endpoint at 
block 1 08, thereby selecting a new re-sampling point as 
a possible clipping endpoint. Execution then loops back 
to block 92. 

[0048] When a point is found that is within the lower 

25 faces at block 92, a similar test is performed at blocks 
110,112, and 1 1 4 to clip the line segment at its opposite 
end. In this case, the "upper* faces of the data volume 
are represented by top face 1 1 6, rear face 1 1 8, and left 
face 120 of data volume 52. It should be noted that the 

30 order in which the segment endpoints are tested is arbi- 
trary, and that the application of the terms "upper" and 
"lower" to refer to different faces of the data volume is 
also somewhat arbitrary, as it depends upon which 
sides of the volume are arbitrarily defined as being "up" 

35 and "down." However, if the line lies entirely outside of 
the data volume, this fact would already have been 
determined in the first clipping test Therefore, there is 
no need to check whether the opposite endpoint of the 
line has been reached in the second clipping test loop, 

40 irrespective of the order in which the endpoints are 
tested. 

[0049] If the point being checked is between the 
lower and upper faces of data volume 52, it is a point at 
which interpolation is performed 122. The interpolation 

45 function is efficiently performed at block 122, because 
no checking need be done in the interpolation function 
to determine whether the point is within data volume 52. 
Next, a check is performed at block 124 to determine 
whether this point is the last point on the segment. If so, 

so execution returns 1 26 to the calling function. Otherwise, 
the next sample point on the line segment is selected at 
block 128, and execution loops back to block 122 until 
the last point is selected. 

[0050] In one embodiment, a feature of this clipping 
55 procedure is that, when all the line segments on a sur- 
face have the same slope, as is the case with those 
aligned with the A axes in Figures 3 and 4, it can easily 
be determined which three faces of volume 52 are 
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required for checking each segment endpoint. Referring 
to a simplified two-dimensional example illustrated in 
Figure 6, a procedure for clipping segment 130 from its 
upper endpoint 132 need only take into account the 
upper Y coordinate face 1 34 of region 1 36 and the lower 5 
X coordinate face 138 of region 136. Similarly, for clip- 
ping segment 130 from its lower endpoint 140, a proce- 
dure need only take into account lower Y coordinate 
face 142 and upper X coordinate face 144. This analysis 
need only be performed once for a surface. In the sim- 10 
plified two-dimensional example, as shown in Figure 6, 
points at which background filling occur are represented 
at 132, 146, 148, 150, 140, 152, 154, and 156. Line seg- 
ment 1 30 is clipped at points 1 58 and 1 60; points 1 62, 
1 64, and 1 66 are representative of points within data 15 
volume 136. In the three-dimensional cases repre- 
sented in Figures 3 and 4, it is recognized that line seg- 
ments 56 and 74 lie entirely outside their respective 
volumes 52. Line segments 60 and 78 are dipped. 
Background filling of output points between points 168 20 
and 170 and between points 172 and 174 of line seg- 
ment 60 in Figure 3 and between points 176 and 178 
and between points 1 80 and 1 82 of line segment 78 in 
Figure 4 occurs. After clipping has been completed for 
any line segment in one embodiment, or for all line seg- 2 s 
ments in another, a specialized sampler is employed. 
[0051] In one embodiment, a specialized sampler 
selected in accordance with a slope of the line seg- 
ments on the re-sampling surface is employed for all line 
segments on the surface. Thus, referring to Figure 3, a 30 
specialized sampler selected in accordance with the 
slope of axis A is employed for that portion of line seg- 
ment 60 on sheet 70, i.e., that portion of line segment 
60 that is between the clipping points 170 "and 174. A 
decision tree is used to select a specialized sampler. An 35 
example of such a decision tree 184 is represented in 
Figure 7. Decision tree 184 is used to select a special- 
ized sampler according to a sign of a slope of line seg- 
ments in a direction of each of the x, y, and z axes. In the 
example shown in Figure 7, an additional choice of a 40 
specialized sampler optimized for each combination of 
possible line segment slopes is provided. Thus, in one 
embodiment, specialized samplers are provided for dif- 
ferent interpolation methods and different sample point 
processing demands. In one embodiment, trilinear inter- as 
polation and nearest neighbor interpolation are used, 
such interpolation methods being known in the art. 
However, the invention is not limited to only these inter- 
polation methods. In one embodiment, other interpola- 
tion methods such as high-order interpolation along one so 
axis of the volume, and lower order interpolation along 
the remaining axes of the volume are provided. Com- 
mon choices of processing demands include inline save 
in interpolated values, or function invocation for each 
sample point for further processing. The second 55 
processing demand is used for planar re-sampling. 
[0052] Figure 8 is a two-dimensional representation 
of an embodiment of planar re-sampling. Extensions of 



this representation to three dimensions will be evident 
to those skilled in the art. In planar re-sampling, a base- 
line segment 186 forms one edge of the sampling 
region 188. From baseline segment 186, sampling line 
segments 1 90 normal to this edge provide scanline seg- 
ments. These scanline segments will generate the inter- 
polated sample values used in an output image. 
[0053] When three different processing or interpola- 
tion methods are provided to implement conventional 
reformat, twenty-four specialized samplers are needed 
for each of the different slope combinations in the deci- 
sion tree example of Figure 7. It may also be desirable 
to support a specialized voxel addressing scheme for 
cases in which the volume dimensions are powers of 2. 
The number of such specialized samplers grows expo- 
nentially with each new variant of re-sampling. There- 
fore, it is desirable to provide a method for automatically 
generating the specialized sampler functions. 
[0054] To provide automatic generation of sampling 
functions, it is desirable to employ encapsulation and 
virtual method techniques in programming the sampling 
functions. This approach provides the simplification and 
reuse benefits provided through the use of object-ori- 
ented paradigms without any degradation of perform- 
ance. Encapsulation is provided, for example, by the 
use of local function variables. 

[0055] Because each specialized sampler is imple- 
mented as one function, local variables are used in one 
embodiment instead of class variables. Various pro- 
gramming languages provide features suitable for gen- 
erating a source code representative of the sample 
functions. 

[0056] In one embodiment written in the C + + pro- 
gramming language, performance is enhanced 
because direct variable references are used instead of 
indirect member variable references. The sampler func- 
tions are generated using class-structured sets of mac- 
ros and a control logic macro. (Class-structured sets of 
macros are those sets of macros having inheritance, 
encapsulation, and virtual methods.) The class-struc- 
tured set of macros define a stepper point class which 
implements the interpolation or processing method. 
This class includes the following methods: 

Const() — This method is the constructor for the 
interpolator/processor. It produces variable decla- 
rations that will be referenced by the other methods. 

lnit(resultBuffer, resultbuffersize) — This initializer 
sets a pointer into the buffer to store interpolated 
values. 

ComputeValue(x, y, z, value) — This method per- 
forms the particular interpolation method for one 
point. 

GetPointNumber(n umber) — This method retrieves 
the re-sampling point number. 0 is the first re-sam- 
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ProcessBackward(value,x,y,z),ProcessFor- 
ward(value,x,y,z) — These methods provide any 
special processing for the value. These methods 
also provide an opportunity to implement the step- 
per for the baseline segment. The forward and 
backward variants allow the processing method to 
track the current point number. 

[0057] In this embodiment, the stepper controller 
macro builds up the overall function by referencing the 
macro-implemented methods of the stepper point class. 
By defining the stepper methods for each style of inter- 
polation or processing, all the specialized sampler func- 
tions are generated. For example, the following code 
segment is used in one embodiment to generate the 
sampler functions using C/C + + program code (the line 
numbers preceding each statement are not part of the 
code but are provided only for reference purposes). In 
another embodiment, equivalent code is written to take 
advantage of corresponding features in another compu- 
ter language. 

1 include(vtkM4ProcessTrilinear.m4) 

2 VtkM4Stepper(RRP,lnline,TR) 

3 VtkM4Stepper(P,P,N,lnline,TR) 

4 VtkM4Stepper(P,N, P, I nline,TR) 

5 VtkM4Stepper(P,N,N,lnline J TR) 

6 VtkM4Stepper(N,P,P,lnline,TR) 

7 VtkM4Stepper(N,P,N,lnline,TR) 

8 VtkM4Stepper(N,N,P,lnline,TR) 

9 VtkM4Stepper(N,N,N,lnline,TR) 
10 

1 1 include(vtkM4Processlnline.m4) 

12 VtkM4Stepper(P,P,Rlnline,NN) 

13 VtkM4Stepper(RRN,lnline,NN) 

14 VtkM4Stepper(RN,Rlnline,NN) 

15 VtkM4Stepper(RN,N,lnline,NN) 

16 VtkM4Stepper(N,P,P,lnline,NN) 

17 VtkM4Stepper(N,RN,lnline,NN) 



19 VtkM4Stepper(N,N,N,lnline,NN) 



w 



20 



21 



22 



23 



24 



15 25 



26 



20 



25 



30 



35 



40 



45 



50 



55 



27 



28 



29 



include(vtkM4ProcessProcess.m4) 
VtkM4Stepper(P,RP,Process,NN) 
VtkM4Steppe r( P, R N , Process, N N) 
VtkM4Stepper(P,N,RProcess,NN) 
VtkM4Steppe r( P, N , N , Process, N N) 
VtkM4Stepper(N,R R Process, NN) 
VtkM4Stepper(N,RN,Process,NN) 
VtkM4Stepper(N,N,RProcess,NN) 
VtkM4Stepper(N,N,N,Process,NN) 



[0058] The first line includes the set of macro defini- 
tions defining the methods for the trilinear interpolation 
style point stepper, The VtkM4Stepper macro is the 
stepper controller macro. This macro uses its argu- 
ments to generate the line-segment slope-specific clip- 
ping code and to create a unique function name. Line 1 1 
redefines the stepper methods with the nearest neigh- 
bor interpolation style point stepper. Line 21 again 
replaces the stepper methods with the process-point 
style point stepper used to implement the baseline seg- 
ment stepper. 

[0059] The point stepper's Const() support good 
run-time performance and encapsulation by generating 
local variable names that are derived from the object 
instance name and the member variable name. Thus, 
data hiding is achieved within any given interpolation 
method. 

[0060] Standard or custom computer systems are 
suitable devices for performing the various method 
embodiments, and constitute inventive apparatus 
embodiments when configured to perform the inventive 
method embodiments. Such computer systems can 
include memory for stored programs and data, one or 
more standard input devices such as an operator con- 
sole, mouse, or joystick, and a video display. A suitable 
computer system is provided as part of the CT imaging 
system in Figure 2. This computer system provides 
mass storage 38, where medical data and a stored pro- 
gram to perform the method may be stored, a computer 
or processor 36, an operator console 40, and a display 
42. Other apparatus embodiments are separate from an 
imaging system, and image data to be processed by 
such embodiments is transferred to the memory of the 
computer in any suitable manner. 
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[0061] Consider a SUN® UltraSPARC-li™ compu- 
ter (Sun Microsystems, Palo Alto, CA) having a 248 
MHz processor. Generation of 512x512 output images 
from a volume of 16-bit voxels with dimensions of 
51 2x51 2x64 can be performed at 26 frames per second 
using nearest neighbor interpolation. Generation of 
512x512 output images from a volume of 16-bit voxels 
with dimensions of 512x512x64 can be performed at 5.5 
frames per second with trilinear interpolation. 
[0062] The above described methods and appara- 
tus provide rapid and efficient re-sampling of a three- 
dimensional data volume to permit interactive viewing 
and visualization of data. THe advantages of the inven- 
tive methods and apparatus are achieved through 
selection of a surface modeled as a set of line segments 
having identical slope, efficient clipping of the line seg- 
ments, and the use and efficient management of highly 
optimized and specialized functions for re-sampling. 
[0063] Although the invention has been described 
and illustrated in detail, it is to be clearly understood that 
the same is intended by way of illustration and example 
only and is not to be taken by way of limitation. In addi- 
tion, the CT system described herein is a "third genera- 
tion" system in which both the x-ray source and detector 
rotate with the gantry. Many other CT systems including 
"fourth generation" systems wherein the detector is a 
full-ring stationary detector and only the x-ray source 
rotates with the gantry, may be used if individual detec- 
tor elements are corrected to provide substantially uni- 
form responses to a given x-ray beam. Moreover, the 
system described herein performs an axial scan, how- 
ever, the invention may be used with a helical scan. 

Claims 

1. A method for rapidly re-sampling a three-dimen- 
sional volume of data (52) to permit efficient inter- 
active viewing and visualization of the data, the 
method comprising: 

selecting a set of parallel line segments (56, 
58, 60, 62, 64, 66, 68) defining a surface (50) 
having a constant slope (dy/dx) in at least one 
direction (A) and intersecting the three-dimen- 
sional volume; and 
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15 



20 



25 



30 



35 



40 



45 



processing the data at points (132, 140, 146, 
150, 152, 154, 156, 158, 160, 162, 164, 166) 
on the parallel line segments to generate a 
reformatted image of the data. 50 



A method in accordance with Claim 1 wherein 
selecting a set of parallel line segments comprises 
the step of clipping ends (168, 172) of the line seg- 
ments so that the clipped line segments (170, 172, 
174) lie wholly within the three-dimensional volume 
(52), and the processing step comprises process- 
ing the data along at points (158, 160, 162, 164, 



55 



166) on the clipped line segments. 

A method in accordance with of Claim 2 wherein 
selecting a set of parallel line segments (56, 58. 60, 
62, 64, 66, 68) comprises the step of selecting a set 
of coplanar parallel line segments. 

A method in accordance with Claim 2 wherein 
selecting a set of parallel line segments (74, 76, 78, 
80, 82, 84, 86) comprises the step of selecting a set 
of parallel line segments in a surface (1 72) having a 
pair of orthogonal axes (A, B), the surface having 
curvature along exactly one (b) of the pair of 
orthogonal axes. 

A method in accordance with Claim 2, 3 or 4, 
wherein clipping ends of the line segments (56, 58, 
60, 62, 64, 66, 68) comprises the steps of: 

determining a first set of faces (94, 96, 98) of 
the volume (52) for clipping a first endpoint 
(172) of each of the line segments and a sec- 
ond set of faces (1 1 6, 1 1 8, 1 20) of the volume 
for clipping a second, opposite endpoint (1 68) 
of each of the line segments; and 

for each line segment, analyzing the first end- 
point of the line segment and the first set of 
faces to clip at a first clipping point (174), and 
analyzing the second endpoint of the line seg- 
ment and the second set of faces to clip at a 
second clipping point (170). 

6. A system (10) for rapidly re-sampling a three- 
dimensional volume (52) of data for efficient inter- 
active viewing and visualization of the data, said 
system configured to: 

select a set of parallel fine segments (56, 58, 
60, 62, 64, 66, 68) defining a surface having a 
constant slope in at least one direction and 
intersecting the three-dimensional volume; and 

process the data at points (132, 140, 146, 150, 
152, 154, 156, 158, 160, 162, 164, 166) on the 
parallel line segments to generate a reformat- 
ted view of the data. 

7. A system (10) in accordance with Claim 6 wherein 
said system is further configured to clip ends (168. 
172) of the parallel line segments (74, 76, 78, 80, 
82, 84, 86) so that the clipped line segments (170, 
172, 174) lie wholly within the three-dimensional 
volume (52), and to process the data at points (1 58, 
160, 162, 164, 166) on the clipped line segments. 

8. A system in accordance with Claim 7 wherein said 
system (10) is further configured to select a set of 
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parallel line segments (74, 76, 78, 80, 82, 84, 86) in 
a surface (72) having a pair of orthogonal axes (A, 
B), the surface having curvature along exactly one 
(B) of the pair of orthogonal axes. 

5 

9. A system in accordance with Claim 7 wherein said 
system is further configured to select a specialized 
sampling function corresponding to a slope (dy/dx) 
of the parallel line segments (56, 58, 60, 62, 64, 66, 

68) and to compute the selected specialized sam- 10 
pling of the data at points (1 58, 1 60, 1 62, 1 64, 1 66) 
along the clipped line segments. 

10. A system in accordance with Claim 7 wherein said 
system is further configured to determine a first set is 
of faces (94, 96, 98) of the volume (52) for clipping 

a first endpoint (172) of each of the line segments 
(56, 58, 60, 62, 64, 66, 68) and a second set of 
faces (116, 118, 120) of the volume for clipping a 
second, opposite endpoint (168) of the line seg- 20 
ments, to analyze, for each line segment, the first 
endpoint of the line segment and the first set of 
faces to clip at a first clipping point (174), and to 
analyze the second endpoint of the line segment 
and the second set of faces to clip at a second clip- 25 
ping point (170). 
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FIG . 2 
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(54) Efficient methods and apparatus for resampling three-dimensional datasets 



(57) In one form, the present invention is a method 
for rapidly re-sampling a three-dimensional volume of 
data (52) to permit efficient interactive viewing and vis- 
ualization of the data is provided in which a set of parallel 
line segments (56, 58, 60, 62, 64, 66, 68) defining a sur- 
face (50) having a constant slope in at least one direc- 
tion (A) and intersecting the three-dimensional volume 
are selected and the data is processed at points (132, 



140, 146, 150, 152, 154, 156, 158, 160, 162, 166) on 
the parallel line segments to generate a reformatted im- 
age of the data. A corresponding apparatus embodi- 
ment is provided. Processing of data can proceed quick- 
ly because it is not necessary to check the parallel line 
segments during processing of the points (158, 160, 
162, 164, 166) on the segments to determine whether 
the points are within a given data volume. 
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